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Abstract 

In complex materials various interactions play important roles in determining the material properties. An¬ 
gle Resolved Photoelectron Spectroscopy (ARPES) has been used to study these processes by resolving the 
complex single particle self energy S(£') and quantifying how quantum interactions modify bare electronic 
states. However, ambiguities in the measurement of the real part of the self energy and an intrinsic inability 
to disentangle various contributions to the imaginary part of the self energy often leave the implications 
of such measurements open to debate. Here we employ a combined theoretical and experimental treat¬ 
ment of femtosecond time-resolved ARPES (tr-ARPES) and show how measuring the population dynamics 
using tr-ARPES can be used to separate electron-boson interactions from electron-electron interactions. 
We demonstrate the analysis of a well-defined electron-boson interaction in the unoccupied spectrum of 
the cuprate Bi 2 Sr 2 CaCu 208 +xCharacterized by an excited population decay time that maps directly to a 
discrete component of the equilibrium self energy not readily isolated by static ARPES experiments. 

PACS numbers: abed 
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INTRODUCTION 


A host of interactions play a role in determining the properties of complex quantum materials. 
Foremost among these are interactions among the electron quasiparticles, and between electron 
quasiparticles and bosonic excitations such as phonons and spin fluctuations. Quantifying these 
interactions is essential to understanding the fundamental and emergent phenomena in complex 
materials. In particular, in strongly correlated materials the phenomenology of these interactions 
remains vital to understanding the behavior of charge and spin density waves, electrical conductiv¬ 
ity anomalies and unconventional superconductivity. However, the analysis is challenging because 
different experimental methods can easily lead to different conclusions. The difficulty primarily 
arises because the interactions overlap in energy, making them difficult to disentangle. 

Interactions quite generally give rise to the quasiparticle lifetime, or equivalently an excitation 
energy linewidth. These lifetimes are described by the energy dependent imaginary part of the 
quasiparticle self-energy, Im T,{E), which represents the effect of the interactions on the quasipar¬ 
ticle at energy E and determines various properties of the material. Each of the interactions has a 
specific contribution to Im S which in equilibrium add to result in a total linewidth according to 
Matthias sen’s rule. This superposition of single particle self energies complicates the disentangle¬ 
ment of the individual interactions, which is desirable to be able to understand the electron-boson 
(e-b) and electron-electron (e-e) interactions. 

Quasiparticle lifetimes or their mass renormalizations can in some cases be directly measured 
using angle-resolved photoemission spectroscopy (ARPES) by studying the linewidth r(E) or 
effective dispersion E{k), obtaining the imaginary or real parts of the self-energy (Im S and 
Re S), respectively [|I]-l5l Recently, femtosecond time-resolved ARPES (tr-ARPES) [l6l-[Tn has 
come to the forefront as another method for studying the quasiparticle lifetimes through analyzing 
population dynamics, and thus ostensibly the interactions. In these experiments, an ultrafast laser 
pump excites quasiparticles from their ground state, and their subsequent energy- and momentum- 
resolved relaxation to equilibrium is probed by analysis of the photoelectron spectrum generated 
by a second, much weaker, time-delayed probe pulse in the UV or XUV spectral range. The 
lifetime analysis of the laser excited population relaxation by kinetic rate equations or density 
matrix formalism jbl Il2l has proven successful for quasiparticle energies between approximately 
0.5 to several eV in a limit where relaxation is dominant at all times and quasiparticle lifetimes 
can indeed be determined in the time domain. iiiiMa 
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At energies elose to the Fermi level E^p the situation is more challenging because two effects 
become important and hamper the lifetime determination. There the various interactions do not 
simply add to determine the population dynamics [fT^ . Here we propose to exploit this seeming 
complication to separate the microscopic interactions of interest using population dynamics ob¬ 
served in tr-ARPES. Immediately after the pump, the primary electronic excitations relax by e-e 
scattering, which leads to a population redistribution of secondary electronic excitations. Subse¬ 
quently, coupling to phonons becomes in general the dominant contribution [|4l|71|. As a result, the 
electron-phonon (e-p) coupling leaves its fingerprint on the population dynamics. 

In this work, we show that microscopic conclusions can be drawn from this particular finger¬ 
print. Using tr-ARPES, we observe population relaxation dynamics which we describe by inter¬ 
action with a specific phonon after excitation by a light field. We apply this analysis to the model 
conducting system Bi 2 Sr 2 CaCu 208 +a;(Bi 2212 ), which is a good candidate because it exhibits pro¬ 
nounced e-p scattering at discrete and well-known energies [|T]-l3l. We find quantitative agreement 
between experiment and theory. By explicit inclusion of the pump light field in our theoretical 
description as the primary electronic excitation mechanism, we avoid the assumption of an ini¬ 
tially thermalized electronic distribution which was required in previous work [|7l. The population 
dynamics are calculated using a numerical time-dependent Green’s function approach. [fTTll The 
photoexcited electrons scatter on an ultrashort time scale by e-e interactions. Because these inter¬ 
actions conserve the total electronic energy in the electron system, they cannot be responsible for 
the dissipation of the excess electronic energy. The dissipation is achieved by electron-phonon (e- 
p) interactions which transfer the energy from the electron system via the hVt 75 meV phonon 
to a heat bath representing acoustic phonons. We validate this assertion by a comparison between 
the tr-ARPES measurements and a theoretical calculation which only includes e-p interactions in 
the relaxation dynamics. 

A measurement of population dynamics on the copper-oxide material Bi2212 in its normal 
state is used both because this material is known to host very strong e-b interactions throughout 
its phase diagram iU [TSl [T3 and to demonstrate the ability of our approach to discern such in¬ 
teractions even when obscured by other physics, as is the case in the normal state of Bi2212. [[T9ll 
The experimentally observed relaxation dynamics exhibit two characteristic timescales, rshort(-£') 
and riong(E), of which Tshort(E) exhibits a sudden speedup at energies above for the lowest 
pump fluence used. The newly discovered short timescale Tshort is found to account for a large part 
of the lifetime obtained from the equilibrium line width tqp{E) measured near Ep along the zone 
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FIG. 1. Time domain separation of electron-electron and electron-boson scattering, a. Illustration 
of the dissimilarity between e-b and e-e scattering for population dynamics. Electron-electron scattering 
(blue arrows) by nature maintains the amount of energy in the electrons, limiting it to rapid redistribution 
of energy within the system. On the other hand, e-b (red arrows) scattering can carry energy out, leaving 
a measurable impression on the population dynamics. This causes the population dynamics to decay with 
the characteristics of the e-b scattering, b. Population decay rate extracted from the time-dependent density 
n{E,t) along the zone diagonal kx = fey as a function of the quasiparticle energy E for various Coulomb 
interaction strengths U. Here, = 0.02 (all interaction strengths are in eV^) and the phonon frequency 
n = 0.1 eV. 


diagonal. By comparison to theory, we show that these time seales ean arise from the eoupling to 
a single phonon, namely the 75 meV O breathing mode[|Tl|^. This diseovery suggests a new use 
for time-resolved measurements —as a seleetive probe of e-b interaetions in materials where their 
presenee and strength are obseured by other proeesses. 

We begin with some general remarks regarding the dynamies of population distributions when 
driven out of equilibrium. The pump deposits energy into the eleetrons. If e-e interaetions are the 
only interaetions present, the exeess energy must remain in the eleetronie system, and although 
the seattering ean eause a rapid ehange in the energy and momentum dependent population, it 
ean only therm ali z e the eleetrons at a higher temperature determined by the absorbed energy. 
However, a eoupling to a bosonie bath of phonons ean draw out the energy by transferring it to the 
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lattice, which has a larger heat capacity than the electronic system. The determining factor for the 
dynamics is then the rate at which e-p scattering dissipates the energy, which depends critically on 
the e-p coupling strength and the phonon frequencies. 

The full dynamics of this process, including e-e and e-p interactions, involves a complex in¬ 
terplay beyond this simple model. Energy relaxation of the photoexcited carrier population is 
expected to occur on three distinct timescales [fT7l (i) electron population redistribution on a fem¬ 
tosecond timescale mediated by e-e scattering; (ii) transfer of energy from the electrons to strongly 
coupled phonons on a tens- to hundreds of femtoseconds timescale; and (iii) full system thermal- 
ization on picosecond or longer timescales, on which the excited population of the strongly cou¬ 
pled phonons decays through anharmonic coupling or transport. In Fig. we illustrate processes 
(i) and (ii), which are the relevant ones for the short term dynamics. The measured and calculated 
time traces typically show a rise at f = 0, followed by a decay. The dynamics are determined by 
a combination of processes (i) and (ii), where e-p and e-e scattering cooperate in a non-additive 
fashion to equilibrate the electrons with each other as well as with the phonon bath. However, only 
e-p interactions can carry energy away from the electrons in a non-magnetic system. 

To show that the e-p scattering is the determining factor in the return to equilibrium, we per¬ 
formed simulations of the pump-probe process using non-equilibrium Keldysh Green’s functions. 
We treat the e-p and e-e interactions self-consistently, with both interactions being evaluated at 
second order in perturbation theory. Here, we go beyond previous work lfTTl |2Tl |23 and include 
the effects of e-e scattering in second order self-consistently. We calculate the time-dependent 
density n(k, t) along the zone diagonal after excitation by a short 1.5 eV pump pulse for a sim¬ 
ple nearest-neighbor tight-binding model and then fit the time traces with a single exponential. 
To illustrate that the resulting decay rates reflect the e-p scattering even when e-e scattering is 
present, we obtain the decay rates for various e-e scattering strengths (U). Fig.[^ shows the decay 
rates obtained for quasiparticles at energy E above Ep. The rates without e-e scattering show a 
step at the phonon frequency ifTTl |22]| . Strikingly, even when the e-e scattering is increased be¬ 
yond the e-p scattering strength (U^ > 0.02), the step remains visible in the data. We also note 
that the high-energy decay rate is smaller when e-e scattering is present, suggesting that there is 
some competition between the two processes. This originates from the fact that the population 
dynamics is a fundamentally different quantity than the dynamics (and therefore lifetime) of a 
singly-excited quasiparticle. The latter is determined by the wavefunction decay, and the former 
by the interactions the quasiparticles with each other and with external phonons. 
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FIG. 2. Equilibrium and excited photoelectron spectra, (a) Photoelectron intensity I{E, A:||, teq) along 
the node of Bi2212 at equilibrium, fgq = -5 ps and equilibrium temperature T = 120 K. b) Difference 
spectrum A/ = I{E, , t) — I{E, , teq) for = 35;uJ-cm“^, t = 50 fs. Light and dark gray bars in panel 
b) depict energy bins AE = 20 meV centered at F^bin used to evaluate /norm(L^bin, t)- 


EXPERIMENT 

We now demonstrate that this analysis ean be applied in real materials by making a quantitative 
eomparison between theoretieal simulations of the type shown above and a tr-ARPES experiment. 
A typical equilibrium spectrum I(E, /cy, teg) along the nodal direction of Bi2212 is shown in Fig. 

Here I is photoelectron intensity and teq = —5 ps precedes the pump pulse. Though the 
dispersion kink is difficult to observe without spectral deconvolution at this resolution ( 23 , the 
increase in coherence (or decrease in linewidth) of the states between —h^l and Ep is still clear. In 
the e-b picture, this is due to a decrease in scattering rate for energies below the boson energy hil. 
Fig. shows the difference spectrum AJ (E, k\\ ,t) = I{E, k\\ ,t) — I{E, fc|| , teq) where f = 50 fs, 

and at incident pump fluence <I) = 35/iJ ■ cm“^. For this fluence, by f = 50 fs electrons excited 
above the equilibrium kink energy +hkl and holes injected below —hkl have largely relaxed, while 
those within the boson window —h^l < E < +h^l persist. Considering that the pump photon 
energy of 1.5 eV (see Methods section for details) is well in excess of this 2hkl range, our tracking 
of the excitations in the vicinity of Ep implies we are primarily observing secondary electrons 
and holes and their relaxation towards Ep, rather than the initial excitation [|T^ . Indeed, Fig. 
shows the largest pump-induced intensity changes within the boson window. The existence of 
this pileup in excited carriers at finite t means the relaxation times for carriers outside the boson 
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FIG. 3. Time-resolved population decay. Intensities /norm(^binj t) at = 315 //J • cm“^ for selected 
binding energies E'bin = 30 meV (blue circles), 70 meV (black circles) and 130 meV (red circles), respec¬ 
tively, on a logarithmic intensity scale. Solid lines are fits to Eq. [T] Data for each E'bin are offset by a factor 
of two for clarity. Dotted lines are guides to the eye extending through /norm (-Ebim t) at times for which 
population decay is dominated by the scattering mechanism responsible for Tshort- Green points are from 
InoTm{E,t) integrated over all /c|| and 1.45 < E < 1.5 eV for t < 200 fs, representing mainly the cross 
correlation (XC) time of pump and probe pulses. The Gaussian fit to the green points (green dashed line) 
yields the XC width. Arrows denote time at which Tiong dominates /norm(f)- 

window must be considerably shorter than inside since the former carriers have already relaxed by 
this time. [l8l|22l Below, we show that the significant change in relaxation times in the unoccupied 
part of the spectrum around +hil results from the sudden change in phase space with increasing 
E for the inelastic decay of excited carriers by emission of a single predominant bosonic mode. 
This conclusion is based upon quantitative agreement between our experiment and the theory of 
fluence-dependent, boson-derived relaxation times which do not account for e-e scattering. 

To quantify these effects, we introduce a procedure to energy-resolve the photoexcited popu¬ 
lation relaxation times r above Ep in a manner that allows comparison to the equilibrium quasi¬ 
particle lifetime tqp determined by analysis of energy- and momentum distribution curves (EDCs 
and MDCs, respectively) in equilibrium ARPES. This is accomplished by dividing the tr-ARPES 
spectrum into discrete energy bins of equal width AE = 20 meV (Eig. |^). Integrating the in¬ 
tensity within bins over k\\ and AE results in time- and energy-dependent intensities I{Eun,t), 
where Ebin is taken at the bin centers. 


8 







- Low <I) -Medium <5 


High <I) - XC 



FIG. 4. Experimental and theoretical fluence-dependent relaxation, a-c) /norm for Sbin = 30, 70 and 

130 meV, respectively for <1> = 35 (blue squares), 105 (black diamonds) and 315 (red circles) • cm“^. 
Solid lines in these panels are fits to Eq.[T] d-f) theoretical /norm for the same Sbin as a-c for field sfrengths 
of 0.05 (blue), 0.15 (black) and 0.5 (red) in unifs of V/a where a is fhe Cu-Cu laffice consfanf. a-f) Green 
dashed lines show fhe Gaussian XC. 

Our subsequent analysis follows from the observation of a pronouneed biexponential deeay of 
photoexeited carrier populations with respect to Eun and $. Typical energy-resolved data for $ = 
315 /rJ ■ cm“^ is shown in Fig. The normalized intensities /norm = /(-Ebin,f)/max[/(Ebin,f)] 
are seen to decay with clearly separable short and long time components (rshort(-Ebin) and 
TongC-E-bin)), the former lasting several hundred femtoseconds. It is this direct observation of 
two exponential time scales within the population decay, rather than one, that allows the present 
analysis of the electron population decay mechanisms in the cuprates Iflbll . Fits to /norm(-Ebin, t) 
at fixed /^binin Fig. as well as all subsequent fits used to determine Tshort and riong, respectively, 
are performed using the function 

/norm(f) = 0(f) (AshortC” + AongC ) (g) R^ (1) 

where 0(f) is the Heaviside function, A’s are intensity amplitudes and is a unit-normalized 
Gaussian of width equal to the cross-correlation (XC) of pump and probe responses at the sample 
surface. 

To facilitate examination of both the energy and fluence dependencies of Tshort and its com¬ 
parison to theory we show representative /norm(f) for several Eun and together with the fits in 
Fig.|^-c. The main trends in the data are the overall decrease in Tshort with increasing Ebin for all 
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fluences and the changes in the ordering of the curves for t > 100 fs, or equivalently the magnitude 

of "Tshort- 

Plotting TshortC-S'bin) for each fluence (Fig. &). we observe a highly unusual relationship be¬ 
tween the fluence dependence of rshort(-Ebin) and the value of F^bin relative to hD.. At the lowest 
fluence rshort(-Ebin) shows a pronounced step at Ml with much shorter times above +Ml 

versus below. As pump fluence is increased, relaxation times below Ml decrease while those above 
Ml increase. The step in Tshort at Ml for 35 /rJ ■ cm“^ is reminiscent of the step in h/rqp (shown 
in Fig.|^)) accompanying boson-induced mass renormalizations observed in ARPES ifTl l4l |2^. 
In ARPES, a single photohole injected above a bosonic mode energy has a shorter lifetime than 
one injected below because an additional interaction channel is open to the higher energy state. 
In energy- and momentum distribution curves (EDC’s and MDC’s, respectively) this phenomenon 
is observed as a step, centered at the mode energy, in the Eorentzian widths of the states. Eike- 
wise, an electron photoexcited above a bosonic mode’s energy will decay faster than one excited to 
empty states below because the additional channel for energy loss due to boson emission is open 
to it.lfTTlI 

The fluence dependence of these phenomena indicate a strong connection between the step 
in Tshort observed above Ep and the nodal kink observed below. It is the process of “filling” 
states within the boson window [fTTII that leads to the characteristic speeding-up and slowing- 
down observed here as well as the weakening of the 70 meV kink structure observed in previous 
experiments ll23l 1^ . As <I) is increased the dynamics of photoexcited carriers change from a 
regime reflecting essentially single-particle physics to one further characterized by the dynamics 
of populations of interacting, hot particles. The parallel evolution, with increasing fluence, of 
the broadening step in Tshort and the previously observed reduction in effective mass of the nodal 
kink below Ep therefore encourages a description of the former effects within the same theoretical 
framework as the latter. The absence of fluence dependence for Tiong(-Ebin) (see Methods) points to 
its likely origin in the coupling of hot electrons to a broad spectrum of collective excitations such 
as acoustic phonon modes. 


Comparison of Experiment to Theory 

With e-p scattering being the dominant mode of energy dissipation, we now focus on a theoret¬ 
ical description of (ii), given experimentally by Tshort. and described by the Holstein model, where 
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FIG. 5. Fingerprint of electron-boson scattering, a) Experimental Tshort(£'bin) for $ = 35 (blue squares), 
105 (green triangles), 165 (black diamonds) and 315 (red circles) /iJ-cm~^, respectively. Dotted blue and 
red lines are guides to the eye for the 35 and 315 data, respectively, b) Theoretical Tshort(f^^bin) 

for field strengths of 0.05 (blue), 0.15 (green), 0.25 (black) and 0.50 (red) and equilibrium tqp = h/2lm S 
(purple). The mode energy at 75 meV is marked by the dashed black lines. Inset) (blue, right axis) Theo¬ 
retical a^F{E) with e-b coupling constant a and dispersion F{E). (red, left axis) The low energy lifetime 
dampening mode at 15 meV is magnified by lOx for clarify. 

electrons are coupled to a single Einstein mode. Since the main feature and overall scale is ex¬ 
pected to be due to the phonon coupling as discussed above, we neglect e-e scattering here. Excited 
carriers relax via their coupling to this reservoir, which is assumed to have infinite heat capacity 
and remains unchanged. The equilibrium tqp{E) in Eig.[^ is defined through the equilibrium line 
width due to a single well-defined phonon mode centered at h^l = 75 meV. The dimensionless e-p 
coupling constant A is chosen to be 0.2 to match the relaxation times of 20 fs at E > hil. 
The same value for A was found to bound the coupling of hot electrons to a particular subset of 
hot phonons in earlier tr-ARPES experiments on Bi2212 [[71, and is similar to that obtained from 
some equilibrium ARPES experiments [^. The model also includes a weakly coupled mode at low 
energies to prevent an infinite tqp within the phonon window {—hQ < E < +}tFI). It should be 
noted that the inclusion of e-e scattering appears to decrease the size of the step in the scattering 
rates, slowing down the decay at energies above hil as can be seen in Eig. This could lead to a 
lower apparent A in time-resolved experiments. 

Representative /norm(-Ebin, f) extracted from theoretically generated time-dependent occupied 
spectral functions are shown in Eig. |^-f for field strengths comparable to the experiment. They 
successfully reproduce the reversal of the fastest rshort(-Ebin) as a function of fluence for Ebin above 
as compared to below +hil. Here, we focus on processes (i) and (ii), and ignore the long-time 
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component, fitting the ealeulated /norm(-Ebin, with a single exponential to produee theoretieal 
'rshort(-£^bin) for several fluenees, Fig.|^. The theoretieal rshort(-Ebin) exhibit the same step at +hVt 
as the experimental data in Fig. |^. The theory largely reproduees the inerease in rshort(-Sbin) for 
Ehin > +hil and deerease in rshort(-Ebin) for 0 < Eun — Ep < +hQ, as a funetion of inereasing flu- 
enee. The remarkable quantitative agreement between the measured and ealeulated /norm(-Ebin, t) 
within the phonon window, Fig. demonstrates that seattering due to phonon emission is suffi- 
eient to deseribe the population dynamies measured with tr-ARPES, allowing for the measurement 
of a single eomponent of the equilibrium self-energy. 


DISCUSSION 

The experiment and theoretieal ealeulations both eapture the dynamies of energy dissipation 
through eleetron-phonon interaetions. The fingerprint of the e-p interaetion through its partieular 
phase-spaee restrietions is seen in the data through both energy and fluenee dependenee of the 
population deeay rates. 

Understanding the detailed nature of the equilibrium self-energy is thought to be vital to under¬ 
standing a host of phenomena in strongly eorrelated materials and beyond, not least those derived 
from non-Fermi liquid properties that would be expeeted to show up most strongly in eleetronie 
ehannels. In equilibrium ARPES, it is diffieult to disentangle phononie from other eontributions 
to the self-energy, espeeially when several different interaetions are present. 

The eombination of eleetron-eleetron (in partieular at higher temperatures), eleetron-boson, 
and eleetron-impurity seattering —just to name a few —eollude and inerease the overall linewidth, 
making it diffieult to separate a single eontribution. With the teehniques demonstrated here, we 
now have the ability to isolate the phonon eontributions direetly from the total equilibrium self en¬ 
ergy through measurement of population dynamies with tr-ARPES with high temporal resolution 
and subsequent eomparison to ARPES data aequired with high energy resolution. 

This finding has important ramifioations beyond the present study. Here, we have foeused on 
a material where the presenee of a eomplex interplay of e-e and e-b interaetions is well known, 
but the general topie reaehes mueh further. This reaeh ineludes other strongly eorrelated mate¬ 
rials sueh as the Ee-based supereonduetors, heavy fermion eompounds, and manganites, as well 
as more weakly eorrelated materials in whieh strong e-p interaetions ean still oeeur. Topologieal 
insulators are a field where the study of population dynamies is eommonly used, and the method- 
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ology presented here ean aid in the interpretation and analysis of the results [|26] - I3^ . Moreover, the 
partieular relevanee of bosonie mode eouplings for exeited state population dynamies whieh is of 
high interest, for instanee, in light-harvesting applieations. We point out that this relevanee is of¬ 
ten implieitly assumed, for example in reeent ab initio eomputations of equilibrium self-energies, 
from whieh information about hot earrier dynamies for solar eells was inferred, ll^ Our study 
provides a firm ground for sueh implicit assumptions, and the methods laid out by us can guide 
further efforts in this direction. 


METHODS 

Sample Preparation and Setup of tr-ARPES Experiment 

Our trAPRES experiment was performed on single crystals of optimally doped Bi2212 (Tc = 
91 K) grown by the floating zone method. Tc was confirmed by SQUID magnetometry. Samples 
were cleaved in situ at the lattice equilibrium temperature, 100 K, in a base vacuum of 5 x 10“^^ 
mbar. Pump pulses of 800 nm wavelength and 40 fs duration, at 250 kHz repetition rate, were 
produced by a Coherent RegA 9040 regenerative Ti:Sapphire amplifier. The 200 nm probe beam 
was produced as the fourth harmonic of part of the RegA’s 800 nm fundamental using nonlin¬ 
ear crystals [13 ■ The pump-probe cross-correlation (XC, Gaussian full width at half maximum) 
was determined to be 100 fs at the sample surface. Photoelectron spectra with an overall energy 
resolution of 55 meV, set by the bandwidth of the laser pulses, were recorded using both angle in¬ 
tegrating time-of-flight (TOE) ||3 and position sensitive TOE (pTOE) [IM1 electron spectrometers. 
The k resolution of the pTOE was 0.0025 and the angular acceptance of the TOE was ±3°. 


Theoretical Modeling of the Data 

Eor the modeling, we use a tight binding parametrization of the one-band model for cuprates 
ll35ll with a filling of 0.42 per spin, i.e. 16% hole doping. The methodology for the calculations is 
a self-consistent Keldysh Green’s function approach as described in Ref. [13 

We couple the electrons in this band to a spectrum of bosons characterized by an Eliashberg 
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function a^F{E) for a spectrum of modes with frequeney VL^, whieh is 
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(2) 


In partieular, we use a dominant boson mode eentered at hVL = 0.075 eV, with = 0.0065eV^ 
eorresponding to dimensionless eoupling A = 0.2, extraeted from the slope of the real part of the 
retarded equilibrium self-energy at zero energy: 


A = 


aReS(E) 

as 


E=Q 


(3) 


To ensure finite relaxation times at low energy, we add a seeond weakly eoupled mode eentered 
at 0.007 eV with g'^ = 2 x 10“®eV^. The frequeney of the strongly eoupled bosonie mode was 
ehosen to mateh the experimentally determined position of the erossing point of fluenee- and 
energy-dependent deeay times, as well as the nodal kink position. The eoupling strength was 
adjusted to give roughly the eorreet time seale in the low-fluenee limit for E > hVt. 

We exeite this metallie single-band model system by a spatially homogeneous pump pulse, 
whose veetor potential eouples to the band eleetrons through the Peierls substitution k —)■ k—eA(f) 
where k is waveveetor, e is the eleetron eharge and A(f) is the time-dependent veetor potential. 
(Here the speed of light c = 1.) The veetor potential A(f), whieh is in the Hamiltonian gauge, 
points along the zone diagonal and has a temporal shape with a Gaussian envelope (cr = 20 fs), 
oseillation frequeney hw = 0.35 eV, and peak strength A^ax- For the simulations in this work, we 
use field strengths = 0.05, 0.15, 0.25, and 0.50 in units of V/oq. For the given parameters 
and a typieal lattiee eonstant Oq = 3.8 A (Cu-Cu distanee), these field strengths eorrespond to 
peak eleetrie fields of 0.65, 2.0, 3.3, and 6.5 MV om~\ respeetively. The time-resolved ARPES 
signal is obtained in a postproeessing step using the formalism deseribed in Refs. [22| and |36] with 
probe pulse of width cxpr = 40 fs. The energy-resolved intensity is obtained by integrating the 
momentum- and energy-resolved data along the Brillouin zone diagonal eut. The deeay times 
shown in the main text (Fig. 4d-f) are extraeted using fits to the normalized intensity ehanges, 
binned in energy windows of 20 meV width like the experimental data and eonvolved with a 
Gaussian envelope funetion of full width XC to render the theoretieal results direetly eomparable 
to the experiment. 

The present modeling and analysis are based mostly on the energy seale (70 — 75 meV), ex¬ 
traeted from the nodal kink position (below Ep) as well as the reported population deeay time 
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step position (above Ep). Importantly, statements about eoupling to bosonie modes with energies 
higher than 100 meV eannot be drawn from our analysis. 


Theoretical comparison of e-e and electron-phonon scattering 

The e-e interaetions are treated self-eonsistently at seeond order with a loeal interaetion strength 
Uo: 


where G^{t,t') is the double-time eontour-ordered Green’s funetion lITTl . Sinee the interaetion 
strength depends on the filling, we define an effeetive U : U = (7gn(l — n) for a better eomparison 
to the e-p eoupling strength. For the eomparison, we used a simplified tight-binding model with 
t = 0.25 eV, t' = 0.075 eV and /x = —0.255 eV. We use a strongly eoupled phonon with = 0.1 
eV and g^=0.02 eV^. It should be noted that the parameters for the eomparison were ehosen 
to make the ealeulations eomputationally feasible and make a elear distinetion between e-e and 
e-p seattering, rather than for a quantitative eomparison to data (as in the other seetions of the 
manuseript). 


Tiong VERSUS Ebin AND 4> 


The long timeseale riong, shown in Fig. was found to earry no fluenee dependenee, within 
error. Interestingly, it was found to nearly merge with, or at least beeome indistinguishable, from 
T’short at high energies, though higher statisties and time resolutions would be required to aseertain 
this with eertainty. Regardless, Tshort is found robustly at all fluenees measured. The fluenee 
independenee of riong plays an important role in our extraetion of Tshort at the lowest fluenee, 
partieularly around i^bin = 70 meV. Beeause at $ = 35/iJ ■ cm“^ the statisties are relatively low 
as Eun inereases mueh past 70 meV and at later times, eonsistent fitting requires riong to be held 
eonstant from higher fluenee data. Relaxing this eonstraint of eourse does not alter the trends for 
'Tshort visible in the raw data, nor does it appreeiably alter the results for the lowest fluenee in e.g. 
Fig. 4a of the main text viz a viz the step around 70 meV. It does introduee a bit more uneertainty 
into the error bars but the position of the step remains. 
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FIG. 6. riong versus Sbin> as defined in the main text, for = 105 (green) 165 (black) and 315 (red) 
/rJ • cm“^. Note the absence of the effect of any energy scale on riong- The vertical dashed line at 75 meV 
denotes the theoretical position of +hQ. 
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